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Abstract. A unique pattern selection in the absolutely unstable regime of driven, 
nonlinear, open-flow systems is reviewed. It has recently been found in numer- 
ical simulations of propagating vortex structures occuring in Taylor-Couette and 
Rayleigh-Benard systems subject to an externally imposed through-flow. Unlike the 
stationary patterns in systems without through-flow the spatiotemporal structures 
of propagating vortices are independent of parameter history, initial conditions, and 
system length. They do, however, depend on the boundary conditions in addition 
to the driving rate and the through-flow rate. Our analysis of the Ginzburg-Landau 
amplitude equation elucidates how the pattern selection can be described by a non- 
linear eigenvalue problem with the frequency being the eigenvalue. Approaching 
the border between absolute and convective instability the eigenvalue problem be- 
comes effectively linear and the selection mechanism approaches that of linear front 
propagation. 

1 Introduction 

In many nonlinear continuous systems dissipative structures branch off a ho- 
mogeneous basic state when the external stress exceeds a critical threshold. 
Examples for these transitions are Taylor-Couette flow, Rayleigh-Benard con- 
vection, binary-fluid convection, flame-front propagation, and some chemical 
or biological processes [1]. Often, for a fixed configuration of parameters and 
boundary conditions (BCs) a continuous or discrete family of patterns with 
different wave numbers is stable. Their stability regime, e. g., a band of wave 
numbers is usually limited by secondary instabilities. The stable structures 
within such a band can be generated by appropriately engineered time his- 
tories of the control parameters and/or by properly changing the BCs, e. g., 
the system size. The most intensively investigated examples in this respect 
are the structures of Taylor vortices [1, 2, 3, 4, 5, 6, 7] in an annulus between 
concentric cylinders of which the inner one rotates, convective roll vortices 
in horizontal layers of one-component fluids [1, 8, 9, 10, 11, 5], and binary 
mixtures heated from below [12, 13]. 

This multiplicity of solutions of the underlying nonlinear partial differen- 
tial equations that stably coexist for a fixed configuration of parameters and 
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BCs seems to disappear in open-flow systems: Recent numerical simulations of 
Rayleigh-Benard convective roll vortices traveling downstream in an imposed 
horizontal Poiseuille flow [14, 15, 16, 17] and of Taylor vortices propagating 
downstream in an axial through-flow [18] showed that the spatiotemporal 
structure of these propagating vortex (PV) flows is uniquely selected, i.e., 
independent of parameter history, intial conditions, and system size in the ab- 
solutely unstable regime. Therein the pattern starting, e. g., from a spatially 
localized perturbation can grow in upstream as well as in downstream direc- 
tion [19]. The structure expansion proceeds until the upstream (downstream) 
moving front encounters in a finite system the inlet (outlet) and adjusts to 
the inlet (outlet) BC. The final pattern resulting in such a situation shows 
a characteristic streamwise profile of the amplitude growing with increasing 
distance from the inlet and of the wave number, and a characteristic global, 
i.e., spatially constant oscillation frequency associated with the downstream 
motion of the pattern. 

By contrast, in the convectively unstable regime [20] initial perturbations 
are blown out of the system — both, the upstream as well as the downstream 
facing front of the growing structure move downstream. Therefore it requires 
a persistent perturbation source like, e. g., noise to sustain a vortex pattern 
in the convectively unstable regime [21]. It should be emphasized that the 
structural dynamics of pattern formation in the convectively unstable regime 
at larger through-flow rates and/or smaller driving rates substantially dif- 
fers from the one reviewed here in the absolutely unstable regime. The latter 
regime is governed by nonlinear contributions in the balance equations; the 
resulting patterns are uniquely selected and insensitive to initial conditions, 
parameter history, and small perturbations. On the other hand, in the convect- 
ively unstable regime the growing patterns are sensitive to initial conditions 
and perturbations. 

In this work we review how such a uniquely selected spatiotemporal pat- 
tern structure in the absolutely unstable regime can be understood as a non- 
linear eigenvalue problem with the oscillation frequency being the eigenvalue 
and the profiles of pattern intensity and wave number determining the cor- 
responding eigenfunction. Although the PV structures selected by the full 
hydrodynamic field equations called Navier-Stokes equations (NSE) for short 
and the Ginzburg-Landau amplitude equation (GLE) approximation differ in 
a characteristic way the selection mechanism looks similar. The frequency se- 
lection seems to result from requiring the spatial variation of pattern envelope 
and phase to be as small as possible under the imposed inlet/outlet BC which 
is analogous to ground state properties of the linear stationary Schrodinger 
equation. 

We also show how the pattern selection process is related to the one oc- 
curring behind a front or "domain wall" that spatially separates an unstable, 
homogeneous state from a stable, structured state. Approaching the border 
between the absolutely and convectively unstable regimes the pattern selection 
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mechanism becomes linear: For driving and through-flow rates on this border 
line the selected frequency is the one resulting from a linear front whose spati- 
otemporal behavior is governed by the fastest growing linear mode. The latter 
is identified by a particular saddle of the complex linear dispersion relation 
over the complex wave number plane [1]. Since the linear dispersion relations 
of NSE and GLE deviate from each other for supercritical control parameters, 
the PV structures selected by NSE or GLE differ in a characteristic way. 



2 The Systems 

We have investigated on the one hand time-dependent, rotationally symmetric 
vortex structures in the Taylor-Couette system with an externally enforced 
axial flow [18] and on the other hand straight convection roll vortices in the 
Rayleigh-Benard system subject to an externally applied plane horizontal 
flow with roll axes perpendicular to the through-flow [14, 15, 16, 17]. In both 
cases a translationally invariant flow state is stable below a threshold of the 
driving, T c (Re) [22, 23, 24, 25] of the Taylor number T or Ra c (Re) [26, 15, 16] 
of the Rayleigh number Ra, respectively, which depend on the through-flow 
Reynolds number Re. 

At these critical thresholds the basic state becomes unstable to PV per- 
turbations via an oscillatory instability and in an infinite system a nonlinear 
PV solution branches off the basic state. We use the relative control parameter 

T 1 Ra 1 m 

a = — 1 or a = — —— - — 1 1 

P T c (Re) P Ra c (Re) V ; 



corresponding to 



T 1 Ra 1 

1 or e = — — — — — 1 (2) 



T c (Re = 0) Ra c (Re = 0) 



to measure the distance from the onset of P V flow for Re ^ and of stationary 
Taylor or convection roll vortex flow for Re = 0, respectively. In this notation 
fi c = 0, i.e., 

£ c^ = T c (Re = 0)- 1 ° r Q (Re) = Ra c (Re = 0) " 1 (3) 

is the critical threshold for onset of PV flow and \i = e/(l + e c (Re)). The shear 
forces associated with the through-flow slightly stabilize the homogeneous 
basic state, so e c (Re) slightly increases with Re [22, 25, 27, 28, 29, 15]. 

In the remainder of this paper we use the Taylor-Couette system as a 
representative example when displaying specific equations. But one should 
keep in mind that our results equally well apply to the problem of propagating 
convection roll vortices in the Rayleigh-Benard setup and also to other similar 
pattern forming systems. 
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2.1 Ginzburg-Landau amplitude equation 

To characterize the secondary PV structure we consider the deviation of the 
velocity field from the basic state as the order-parameter field. Close to the 
bifurcation threshold of PV flow, i.e., for small /i, the flow has the form of a 
harmonic wave. For example, the axial component of the vortex field in the 
Taylor-Couette system reads 

w (r, z; t) = A {z, t) e'^-^w (r) + c.c. (4) 

with a complex amplitude A(z,t) that is slowly varying in the axial direction 
z and time t. Here and in the following all quantities are properly reduced 
[18, 14] to be non dimensional. The critical wave number k c , frequency uj c 
[23, 24, 22, 25, 27, 28], and eigenfunction w(r) [23, 30] depending on the 
radial coordinate r have been obtained from a linear stability analysis of the 
basic flow state as functions of Re. The complex vortex amplitude A(z,t) is 
given by the solution of the ID complex GLE, 

r (i + v g A'^ = a* (1 + ic ) A + ^ (1 + id) A" - 7 (1 + »c 2 ) \A\ 2 A. (5) 

Dot and primes denote temporal and spatial derivatives in the z-coordinate, 
respectively. As a consequence of the system's invariance under the combined 
symmetry operation {z — » —z, Re — » —Re} the coefficients To, £q, 7 are even 
in Re while the group velocity v g and the imaginary parts Co, c\, c 2 are odd in 
Re [16, 22]. An analogous expression of (4) and GLE also holds for convective 
roll vortices in the Rayleigh-Benard system. 

It should be noted that the control parameter range of \i over which (4) 
gives an accurate description, say, on a percent level is indeed very small in 
the Taylor-Couette system: The asymmetry between radial in- and outflow 
intensities rapidly grows with \i and causes higher axial Fourier contributions 
~ e mkz [31, 32, 33] to the velocity field that are discarded in the fi — > 
asymptotics of the GLE approximation (4). In the Rayleigh-Benard system 
the situation is more favorable for a GLE approximation. In any case however, 
the modulus of the first Fourier mode of the vortex structures agrees in both 
systems for Re = as well as for Re 7^ quite well with the one predicted by 
the GLE. 

2.2 Absolute and convective instability 

For small e and Re the control parameter plane is divided into three stability 
regimes - cf. Fig. 1 - characterized by different growth behavior of linear per- 
turbations of the basic flow state. In the presence of through-flow one has to 
distinguish [19] between the spatiotemporal growth behavior of spatially local- 
ized perturbations and of spatially extended ones. Below the critical threshold 
e c (Re) for onset of PV flow (dashed line in Fig. 1) any perturbation, spatially 
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Fig. 1. Stability domains of the basic state. Numerical simulations of the 
Taylor-Couette system have been performed for the parameters marked by sym- 
bols. Dashed line is the critical threshold for onset of extended PV flow and full line 
the boundary (6) between absolute and convective instability. 



localized as well as extended, decays. This is the parameter regime of absolute 
stability of the basic state. 

For e > e c extended perturbations can grow. A spatially localized per- 
turbation, i.e., a wave packet of plane wave perturbations is advected in the 
so-called convectively unstable parameter regime faster downstream than it 
grows — while growing in the comoving frame it moves out of the system 
[19, 20, 21, 34]. Thus, the downstream as well as the upstream facing intensity 
front of the vortex packet move into the same direction, namely, downstream. 
A spatially localized perturbation is blown out of any system of finite length 
and the basic state is reestablished. It therefore requires a persistent perturb- 
ation source like, e. g., noise to sustain a vortex pattern in the convectively 
unstable regime [21]. 

In the absolutely unstable regime (shaded region in Fig. 1) a localized 
perturbation grows not only in the downstream direction but it grows and 
spatially expands also in the upstream direction until the upstream propagat- 
ing front encounters the inlet in a finite system. The final pattern resulting in 
such a situation shows in the downstream direction a characteristic intensity 
profile under which the PV flow develops with increasing distance from the 
inlet. Approaching the boundary between absolute and convective instability 
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the growth length of the PV structure diverges in the absence of any perturb- 
ation source and the PV pattern is blown out of the system. However, in the 
presence of noise there is a transition to a noise sustained structure with a 
characteristic finite growth length depending on noise properties and control 
parameters. 

Within the framework of the amplitude equation the boundary (full line 
in Fig. 1) between absolute and convective instability is given by [21] 

c T V g 

^ conv = 4^ 2 (l + c?) (6) 

corresponding to t% onv = e c + (1 + c^Hconv Thus the absolutely unstable 
regime (shaded region in Fig. 1) is characterized by /j, > H% onv or, equivalently, 
by the reduced group velocity 



wa + c?)^ V A* 



t t u o r'conv /<-7\ 

V g = Z . v 9 = 2 \l— ( 7 ) 



being smaller than 2. It should be mentioned that the GLE approximation (6) 
of iJ. c conv describes the boundary between absolute and convective instability 
resulting from the NSE [28, 15, 16, 17, 35] very well for the small Reynolds 
numbers considered here. 

When comparing results obtained for different /j, and Re from numerical 
simulations of the NSE with those following from the GLE we found it some- 
times advantageous to present them as functions of the reduced group velocity 
V g (7). The deviation 

2 - V g = 2 - ^ii% on Jii (8) 

is a scaled distance from the boundary between the absolutely and convect- 
ively unstable regimes. 



3 Propagating Vortex Patterns 

We have performed numerical simulations of the full 2D NSE for rotationally 
symmetric Taylor vortices and straight convective roll vortices in systems of 
finite length [18, 14, 15, 16, 17]. They revealed in the absolute unstable regime 
for various end conditions patterns of vortices propagating downstream under 
a stationary intensity envelope after transients have died out. The oscillation 
frequency u> of the PV flows is spatially constant while the local wave number 
k, the phase velocity v p = u>/k, and the vortex flow intensity varies. The 
frequency and the spatial variation of the PV pattern depends on the control 
parameters and on the BCs but not on parameter history or initial conditions. 

For simplicity we here only discuss vortex suppressing BCs that are real- 
ized by imposing the homogeneous basic state at the inlet (z = 0) and outlet 
(z = r = 50) of the system. In the GLE approximation this amounts to 
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10 20 30 40 50 

distance from the inlet 

Fig. 2. Spatiotemporal structure of PV flow. Thin lines show vertically displaced 
snapshots of the vortex velocity field at successive, equidistantly spaced times. Thick 
lines show the stationary envelope. The BC at the inlet (left) and outlet (right) 
suppresses any vortex flow there. 

requiring A = at the inlet and outlet. To illustrate the global properties 
of the flow patterns we present in Fig. 2 a hidden-line plot of the axial ve- 
locity field w in the Taylor-Couette system. Thin lines show snapshots of w 
at the radial position r = r\ + 0.225 obtained from the NSE at successive, 
equidistantly spaced times. After transients have died out the vortices propag- 
ate downstream under a stationary intensity envelope (thick line) determined 
by the temporal extrema of w(r\ + 0.225; t) at any z-position. The intensity 
variation of PV flow under the fronts causes there a spatial variation in the 
local wavelength A(z) and in the phase velocity v p (z) = ujX(z)/2ir with the 
oscillation frequency u> of PV flow being globally constant [18, 14]. 

For further characterization of the PV flow structure we found a temporal 
Fourier decomposition of the time-periodic fields, e.g., 

w(r,z;t) = J2^n(r,z)e-^ t (9) 

n 

to be useful. Note that the long-time solution of the GLE 



A(z,t) = iJ(z)e*W*)- m ] 



(10) 
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oscillates harmonically with frequency Q under a stationary envelope 

R(z) = \A(z,t)\ (11) 

after transients have died out. Therefore, the GLE fields contains no temporal 
Fourier mode other than n = 1 whereas the solution of the full NSE fields 
have higher harmonics and a stationary contribution of a zeroth fourier mode 
[18, 14]. 

Sufficiently away from the inlet and outlet we observe a bulk region of non- 
linear saturated PV flow with spatially uniform amplitude and wavelength. In 
the bulk region the temporal modes obtained for finite through-flow from the 
NSE grow for small /j, proportional to / u n/ ' 2 with relative corrections propor- 
tional to [i. Thus, they show the same growth behavior with /j, that the spatial 
Fourier modes [31, 36, 37, 32, 33] of stationary vortices without through-flow 
show as a function of e. The reason is that all fields in the PV state have the 
form of propagating waves 

f b (r,z;t) = f b (r,z-^t) (12) 
in the bulk region - denoted by a subscript b. 




Fig. 3. Scaled growth length L (13) of PV structures versus scaled group velocity 
V g (7). Symbols represent lengths obtained from the NSE for different combinations 
(cf. Fig. 1) of Re and e = 0.0288 (circles), e = 0.114 (diamonds), and e = 0.186 
(triangles). The line shows the scaling behavior of the GLE. 
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With increasing through-flow this bulk pattern is pushed further and fur- 
ther downstream. The growth length / from the inlet over which the amplitude 
reaches half its saturation value depends on /j, and Re. It increases and finally 
diverges when the control parameters /i, Re approach the absolute-convective 
instability border /j™"*, that is when V g (7) reaches 2. In Fig. 3 we compare 
the scaled growth length 



computed from the amplitude equation (solid line) with results from the NSE 
(open symbols) for various combinations of e and Re. Due to the scaling 
property of the GLE [15, 16, 17] keeping in mind the smallness of the imaginary 
parts Cj all values for L obtained by the GLE fall onto one curve for different 
control parameters. The NSE results lie very close to the GLE curve. This 
also holds for propagating convective roll vortices [14, 15, 16, 17]. 

4 Pattern selection 

Our investigations of the PV structures in the Taylor-Couette system [18] and 
in the Rayleigh-Benard system [14, 15, 16, 17] have shown that the through- 
flow causes a unique pattern selection such that the selected structure is in- 
dependent of history and initial conditions and depends only on the control 
parameters and BCs. Thus, the flow causes possible wave numbers within 
the Eckhaus-stable band of stationary Taylor or convection vortex patterns 
without through-flow to collapse to only one uniquely selected PV structure 
for nonvanishing through-flow rates. 

4.1 Pattern selection within the GLE 

Here we elucidate how the pattern selection mechanism of the GLE can be 
understood as a nonlinear eigenvalue/boundary- value problem where the fre- 
quency of the PV pattern is the eigenvalue. Thus, the selected PV structure 
is characterized and determined by the combination of eigenvalue and corres- 
ponding eigenfunction. 

The eigenvalue problem We look for solutions of the GLE of the form 



(13) 




(14) 



q(z) = k(z) - k c = <p'(z), 



(15) 



and constant frequency 



Q = u> — u>, 



(16) 
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We have numerically solved the GLE and verified that its long-time solution 
is indeed of the form (10). Inserting this solution ansatz into the GLE one 
obtains the nonlinear eigenvalue problem 

r (-if2a + v g a') = /j (1 + ic ) a + (1 + ici) a" -7(1 + ic 2 ) \a\ 2 a (17) 

or, equivalently, 

iT (-O + v g q)R + T v g R' = \i (1 + ic )R + 

^(l + i Cl )(R" - q 2 R + iq'R+2iqR l )- 1 {l + ic 2 )R 3 (18) 

as a solvability condition with Q and a(z) being the eigenvalue and associated 
eigenfunction, respectively. Here we consider the BCs Ri n ,out = at the inlet 
and outlet. These requirements fix four BCs — namely, JL4j' notlt = QAi n out = 
— which are necessary to solve the GLE. 

We are interested in solutions for which the variation of R(z) and in par- 
ticular of q(z) is "as small as possible". Such a solution type seems to be 
connected with the eigenvalue Q that is closest to zero. While we cannot 
make statements about the eigenvalue spectrum our numerical solutions of 
the GLE indicate that initial conditions always evolve into a pattern with 
smooth amplitude a(z) and small frequency Q. In that sense the eigenvalue 
problem (17) resembles ground state problems of the Schrodinger equation 
— increasing the spatial variation of the wave function (amplitude a) implies 
higher kinetic energy (frequency f2). 



Inlet behavior While the eigenfrequency Q is in general determined by 
global properties and not by the local variation of the eigenfunctions R(z) and 
q(z), say, at the inlet or in the bulk, it is informative and useful for our further 
discussion to list below relations between Q and structural properties at the 
inlet and in the bulk of the system. From (18) one finds that R = at the 
inlet and outlet implies the relation 

Rin,out = [(2(1 1 • 7 _ ^lin ,out\Ri n ,out 
Co I 1 T * C 1 ) 

at z = 0, r . Hence, whenever the modulus grows with finite slope, R' in out ^ 0, 
the BC R = fixes the wave number at the inlet and outlet to the value 

e _ 1 T ° v a - c i /— c / on >, 

C,0lin,out — — TTCl , 2T— , Vacant; ( ZU ) 

I 4o(l + c 1 ) Vl + cf 

which depends only on Re. This value ensures that the imaginary part of the 
expression in the square brackets in (19) is zero. For Re = one obtains 
qi n out = as Cross et al. [38]. In addition the condition R = yields via the 
real part of (19) two relations 



f>" — T ° V 3 pj 

^in.out C2(i 1 2\ in, out V / 

Col 1 T C l) 
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between the different slopes R' in ^ R' ou t alL< ^ curvatures R'/ n ^ R'out °f the 
modulus at the two respective boundaries that also hold independently of pt. 

Furthermore, one can derive a relation between the eigenfrequency f2, the 
wave number qi n ,out, and its slope q' in out at the inlet and outlet. To that end 
we consider the spatial derivative of (18) at the inlet and outlet with R = 0: 

R" 

ir {Q - v g q) + fi(l + ic ) - £q(1 + ici)q 2 - [t v 9 - 2ig^(l + ic i)]^T 

T3lH 

= Zl{l + ic 1 ){— + Siq'). (22) 

This allows one to solve for -^77- and q' separately. Using (20) and (21) one 
obtains the relation 

r i? = (ci - c )// + T v g q in)OUt - 3(1 + c\)ilq' in out . (23) 

Note that this equation implies q' in = q' out . 

The relation (23) also demonstrates the nonlocality of the pattern selection 
mechanism and of the eigenvalue problem. At first sight one might be tempted 
to infer from the above relation between Q and g 8 -„ , q' in that the latter two fix 
the frequency. However, Eq. (23) holds equally well for q ou t, q'out an d the outlet 
properties do not fix Q either. So the correct interpretation of Eq. (23) is that 
the eigenvalue f2, being the characteristic global signature of the pattern, fixes 
the local quantities q' in out for the given BC. 



Bulk behavior The wave number q(z) away from the inlet and outlet in 
general differs from qi n (= q ut)- In order to obtain the full spatial profiles of 
the eigenfunctions R{z) and q(z) belonging to the eigenvalue Q one has to 
solve the eigenvalue problem (17) with the BC Ri n ,out = 0. Alternatively one 
can solve numerically the time-dependent GLE. 

When the system size and the control parameters are such that the PV 
pattern forms a bulk part with a homogeneous modulus R and wave number 
q , i. e., where 

R' b = < = q' b = 0, (24) 

then the dispersion relation 

r i? = (c 2 - Co)// + T v g q b + (ci - c 2 )^o?6 (25) 

provides a relation between the frequency eigenvalue f2, the bulk wave number 
q b , and the bulk modulus 

Rl = (26) 

7 

The equation (25) establishes together with (23) also a relation between q b on 
the one hand and g 8 -„ = q out and q' in = q' out on the other hand. 
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Approaching the convective instability boundary (6) we found that q' in 
decreases to zero so that the frequency eigenvalue becomes 

T 0^conv = ~( c + Cl)t* C conv ( 27 ) 

Then the bulk wave number approaches according to (27) and (25) the limiting 
value 

io{lb) C conv — V^conv (28) 

C\~ C'2 

while the inlet wave number rji n is given by (20). Due to the very small 
imaginary parts c 8 - the solution (28) with the plus sign has to be discarded, 
since it corresponds to unphysically large wave numbers. 



Comparison with front propagation It is suggestive to compare proper- 
ties of PV patterns like the one in the Fig. 2 with those behind an upstream 
facing front that connects in an infinite system the basic state at z = — oo 
with the developed PV state at z = +oo. A linear growth analysis of the GLE 
along the lines of [Sec. VI B 3;1] shows that the far tail of a linear front of 
PV perturbations moves with the front velocity 



v F = v g -2^\L{l + cl). (29) 



In the convectively unstable regime < [i < ^% onv the front propagates down- 
stream so that vp > 0. In the absolutely unstable regime /j, > [i c conv it moves 
upstream, i.e., vp < 0. And at the boundary /j, = H% onv the front is stationary 
— vp = 0. We compare the properties of such a stationary front at /j, = H% onv 
with those of the solution A(z,t) (10) with stationary modulus in a semi- 
infinite system (0 < z < oo). The reason for restricting the comparison to 
l^conv i s that an upstream moving front in the absolutely unstable regime will 
be pushed against the inlet whence the spatiotemporal structure of the free 
front gets modified by the inlet. 

In the linear part of the stationary front the flow amplitude 

A(z,t) ~ e *'K^-C(Q s )t] (30) 

varies with a wave number qp = 5RQ S , a spatial decay rate 9Q S , a temporal 
growth rate QC(Qs) = 0, and an oscillation frequency Qp = 3^C(Q« ) • Here 

1 i t i + ci 



CoQs = -2TTi^To Vs = '^/T+W^ (31) 

is the saddle position of the complex dispersion relation 

toC(<3, AO = T v g Q + i(l + ic )n - i(l + ici)£, 2 Q 2 (32) 
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of the linear GLE in the complex Q-plane determined by the condition 

d((Q) 



conv ' 



Q s " ° (33) 



for vp = 0, i.e., at /j, = \f c 

Thus, the stationary front resulting from the linear GLE is characterized 
by the local wave number in the far tail 

S,oqF{v F = 0) = Cl ^li c conv (34) 

V 1 + C i 

and the oscillation frequency 

t (2 f {vf = 0) = -(c + ci)n c conv . (35) 

So the inlet wave number rji n (20) agrees with qp (34) and f2 c conv (27) coincides 
with Qp (35). Thus, right at the border H% onv of the absolutely unstable 
regime the frequency eigenvalue of the PV pattern growing in the downstream 
direction from the inlet at z = agrees with the oscillation frequency selected 
by the stationary, linear front. But the PV front pattern develops from the 
basic state with an exponential growth rate k of the modulus. This leads for 
R — > to R' = kR,R" = k 2 R, ... — > for z — > — oo. On the other hand, 
the PV flow intensity in the semi-infinite system in general drops to zero 
at the inlet with finite R' and R" . Only at the border line /j, = H% onv the 
difference disappears since there R' , R" , ... — > for z = 0. In the absolutely 
unstable regime, V g < 2, the frequency eigenvalues Q = u> — uj c of the GLE 
lie slightly above the limiting value f2 c conv at the border of the absolutely 
unstable regime. The half-tone symbols in Fig. 4 show for three different e how 
Q obtained numerically by integrating the GLE approaches with increasing 
Re the front frequency 3^C(Q« ) = &f(vf = 0) = Q c conv (dashed line). For 
convenience, the merging points with this limit line are marked for the three 
e- values investigated here by small half-tone symbols. 

Inserting at ^ c corlv the linear front frequency Qp(yp = 0) = Q c conv into the 
bulk dispersion relation (25) leads to an expression for the bulk wave number 
far behind the front. For vp = one obtains 



£o(qb)F = Co(qb) C C onv = V^conv (36) 

C\~ C'2 

As an aside we mention that the nonlinear front solution of Nozaki and Bekki 
[Eqs. 2-7;39] yields for our GLE a front that is not stationary at [i c conv but 
rather moves with velocity 




vnb = ~(3W , n 2 ~ l ) v g ( 37 ) 
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Fig. 4. Selected oscillation frequency of PV flow in the Taylor-Couette system ob- 
tained from the NSE (open symbols) and GLE (half-tone symbols) versus Reynolds 
number. The full (dashed) line is the front frequency Qf = ^C(Qs) of the linearized 
NSE (GLE) at the border between absolute and convective instability reduced by 
the viscous diffusion rate is/d 2 across the gap. Parameters are e = 0.0288 (circles), 
e = 0.114 (diamonds), and e = 0.186 (triangles). The small symbols marking Qf 
for these three e values show that the eigenfrequencies of the nonlinear equations in 
the absolutely unstable regime approach the Qf limit curve at the right places. 



in the upstream direction with a bulk wave number 

(c 2 - ci)V8 + 9cf 
[3(1 + c\) + sign{c 2 - ci) x ^/8{c 2 - Cl ) 2 + 9(1 + Cl c 2 ) 2 ] . (38) 

So we conclude that the nonlinear front solution of Nozaki and Bekki is unre- 
lated to the PV patterns at [i c conv . Furthermore, it was noted [40] that localized 
initial perturbations did not evolve into the nonlinear front solution of Nozaki 
and Bekki [39, 41, 42]. However, v^b and {qb)NB differ for small through-flow 
only slightly from the respective values vp = and (qb)F- For example, at 
Re = 2 one obtains v^p = — 0.06i> s and (qb)NB = 1-15(<7{>).F- 

Scaling properties Scaling length and time according to 

z ~ t 

z = v^t ; 1 = ( 39 ) 

CO TO 
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the GLE for the reduced amplitude A = A\J ^ / ^ no longer contains 7 and [i 
explicitly but the reduced group velocity V g (7) and the coefficients Co, c\, C2. 
Thus the reduced selected frequency Q = / '/i and the selected bulk wave 
number % = CoQb/y/]^ depend not only on V g but via co,ci,C2 also on the 
Reynolds number. This i?e-dependence can alternatively be seen — via the 
i?e-dependence of V g and of the c 8 -'s entering (7) — also as an additional \i- 
dependence. The latter dependence is sufficiently strong to prevent a scaling of 
Q with V g alone. However, by reducing the eigenfrequencies C2(e, Re) selected 
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Fig. 5. Frequency shift to — lo c of PV flow in the Taylor-Couette system versus scaled 
group velocity V g (7). Both, NSE (open symbols) and GLE (half-tone symbols) 
results are scaled by the corresponding GLE frequency fl c conv (27) at the convective 
instability border for the respective parameters e = 0.0288 (circles), e = 0.114 
(diamonds), and e = 0.186 (triangles). The small symbols have the same meaning 
as explained in Fig. 4. 



for our three different e-values by the limiting values f2(e, Re c conv ) = fi^onv 
(27) at the border of the absolutely unstable regime we effectively have elim- 
inated the //-dependence and all GLE data (half-tone symbols in Fig. 5) fall 
onto one curve. A similar scaling holds for the bulk wave number 55 divided 
by ( g6 )< onu (28). 
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4.2 Pattern selection within the NSE 

Remarkably enough, by reducing the NSE frequencies with the GLE frequency 
f^conv (27) the NSE results (open symbols in Fig. 5) almost show this one- 
variable scaling with V g . However, for the GLE (u> — u> c ) / f2 c conv increases 
monotonously with the scaled group velocity V g whereas the NSE results show 
a monotonous decrease with V g . This discrepancy between NSE and GLE 
results is caused by the different dispersion relations of the full hydrodynamic 
equations and the approximate GLE — cf. below. 

The NSE wavelengths that are selected when the basic state is enforced 
at the inlet and outlet weakly decrease with increasing Re and e. For this 
BC the wavelengths are mostly smaller than the critical one with a deviation 
of up to 2% for the largest e = 0.186 that was investigated [18]. Qualitative 
and quatitative similar results were obtained for transversal Rayleigh-Benard 
convection rolls in horizontal shear flow [15, 16]. These deviations are much 
stronger than the decrease of the critical wavelength with increasing Reyn- 
olds numbers. To summarize: the wavelengths resulting from the NSE differ 
systematically from those resulting from the GLE. The former decrease with 
growing Re and e while the latter increase. Also experiments [43, 44] yield 
wave numbers which differ in a to our NSE results common distinctive way 
from the proper GLE result (36). Nevertheless we think that the selection 
mechanism is quite similar — cf. the discussion below. 

4.3 Comparison with front propagation — NSE versus GLE 

Let us first consider Taylor vortices without through-flow. Then the frequency 
Q is zero, the GLE contains no imaginary coefficients c 8 -, and the first-order 
spatial derivative is absent since v g = 0. In this case the GLE yields for 
propagating fronts [45] as well as for stationary patterns in finite systems 
[38] the critical wave number, q(z) = 0, all over the extension of the pattern. 
The BC A = enforces the collaps of the supercritical band of stable bulk 
wave numbers of nonlinear vortex patterns to q = 0. This nonlmearly selected 
wave number q = happens to be for Re = the same as the wave number of 
maximal linear growth under a linear front whose spatiotemporal evolution 
is governed by the dispersion Cgle(Q, -Re = 0) (25) of the linear GLE. In 
the presence of through-flow, however, the wave number of maximal growth 
under the linear front — e.g. qp (34) — of the GLE differs from the one in 
the nonlinear bulk far behind the front [(<?{>) f (36)]. 

Now, already without through-flow the dispersion relation of the linear 
NSE differs significantly from that of the GLE so that the wave number of 
largest temporal growth under a linear NSE front deviates from the GLE 
result q = 0. Moreover it has been observed experimentally [3, 10] and nu- 
merically for Taylor vortices [6, 46] or Rayleigh-Benard convection rolls [11] 
that the wave number of real fronts is spatially varying in contradiction to the 
GLE picture. 
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Now let us consider the through-flow case with a nonvanishing frequency 
eigenvalue Q. Then also in the GLE approximation the selected flow patterns 
reveal spatially dependent wave number profiles connecting the wave number 
of maximal growth in the far tail of the front or at the inlet to the one in 
the bulk. In Fig. 4 we compare the eigenfrequencies of the nonlinear NSE 
(open symbols) with those (solid line) that would be selected by a linear front 
that is stationary at the border of the absolutely unstable regime. To that end 
Recktenwald and Dressier [35] have determined the dispersion relation C(Q) 
of linear perturbations ~ e l ^ kc+l ^' ,z ~^ Wc+ '= ' ,t ^ in the complex Q-plane resulting 
from the NSE with a shooting method as described in Ref. [28]. The solid line 
in Fig. 4 represents the frequency fip = mC,(Q s ) selected by the stationary 
front (vp = 0) of the linear NSE. Again, as discussed for the GLE in Sec. 4.1, 
this frequency is determined by the saddle Q s of C(Q) f° r which the temporal 
growth rate QC(Qs) = 0. This defines the boundary between absolute and 
convective instability. With increasing Re the eigenfrequency of the nonlinear 
NSE approaches the frequency selected by the linear stationary front at the 
border of the absolutely unstable regime as indicated for the three e- values 
shown in Fig. 4. The front frequencies for these e- values are marked by small 
open symbols on the curve 3^C(Q« ) t° indicate that the eigenfrequencies of 
the NSE in the absolutely unstable regime (large open symbols in Fig. 4) do 
indeed end at the right positions on this curve. 

This NSE behavior is very similar to the GLE behavior (half-tone symbols 
and dashed line). Note, however, that the linear dispersion relations of NSE 
and GLE are different: the front frequency ^Cnse (Qs) (solid line in Fig. 4) 
is positive and increases with through-flow while %ICgle(Q s) (dashed line) 
is negative and decreases with through-flow. This difference in the disper- 
sion relations seems to be the major cause for the differences in the selected 
patterns. 

5 Conclusion 

We have reviewed propagating vortex structures in systems of finite length 
that propagate downstream in an externally applied flow. Within the abso- 
lutely unstable parameter regime a unique pattern selection is observed. The 
selected PV patterns are independent of parameter history, initial conditions, 
and system length provided the latter is large enough to allow for a saturated 
bulk region with homogeneous pattern amplitude and wave number. But they 
depend on the inlet and outlet BC. The outlet BC influences the structures 
only locally in the vicinity of the outlet while the globally constant frequency 
and the streamwise profiles of pattern amplitude and wave number and their 
bulk values are basically determined by the inlet BC. The selected frequency 
u> of the pattern oscillation is an eigenvalue of a nonlinear eigenvalue problem 
with the set of time periodic fields being the associated eigenfunction. The 
analysis of the appropriate GLE approximation shows that the eigenfunction 
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associated to the selected eigenfrequency varies as smoothly and as little as 
possible under the imposed BC. This property is similar to the ground state 
behavior of the linear stationary quantum mechanical Schrodinger equation. 

Different lateral BCs entail different frequencies and eigenfunctions, i.e., 
lateral pattern profiles. For example, inlet conditions that exert a phase pin- 
ning force on one or more of the fields of traveling vortices reduce the fre- 
quency. The spatial growth behavior and the pattern profile of the PV struc- 
tures in the region between inlet and bulk saturated pattern is strongly influ- 
enced by the inlet BC. But with increasing growth length of the pattern from 
the inlet the influence of the latter on spatiotemporal properties of the convec- 
tion pattern becomes weaker. While the frequencies and bulk wave numbers 
obtained for various BCs and through-flow rates well inside the absolutely un- 
stable regime differ substantially they all fall onto a linear dispersion curve. 
Finally, when approaching with increasing Re the border between absolute 
and convective instability the eigenfrequencies and also the spatial pattern 
profiles belonging to different BCs approach each other. On this border the 
eigenvalue problem becomes effectively linear and the pattern selection mech- 
anism becomes that one of linear front propagation. 
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